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Abstract — A simple decoding procedure for algebraic-geo- 
metric codes CulD.G) is presented. This decoding procedure 
b a generalization of Peterson's decoding procedure for the 
BCH codes. It can be used to correct any [(il' - 1 )/2J or fewer 
errors with complexity 0(n ' J, where </* is the designed minimum 
distance of the algebraic-geometric code and // is the codelength. 
Index Terms— Error-correcting codes, algebraic-geometric 
I procedure, correcting [(./' - 1)/2J errors. 



I. Introduction 

THE MOST important development in the theory of error- 
correcting codes in recent years is the introduction of 
methods from algebraic geometry to construct linear codes. 
These so called algebraic- geometric codes were introduced 
byGoppa. In 1982, Tsfasman, VISdut, and Zink [1] obtained 
an extremely exciting result: the existence of a sequence of 
codes that exceeds the Giibert-Varshamov bound [2], For 
this paper, they received the IEEE Information Theory Group 
Paper Award for 1983. Since then, many papers dealing with 
algebraic-geometric codes have followed |3]-[10]. 

Gtx)d code constructions are very important. Moreover, it is 
desirable and important to derive simple decoding procedures 
which can correct as many errors as possible. Justesen el 
til J 1 1 ] first presented a decoding procedure for codes from 
nonsingular plane algebraic curves. This decoding procedure 
can only correct {{<T -g- 1)/2J or fewer errors, where d* 
is the designed minimum distance of the code and <j is the 
genus of the curve involved in the construction. Skorobogatov 
and VladuJ, [ 1 2] generalized their ideas and gave a decoding 
procedure which can correct any [(rf* - g - 1)/2J or fewer 
errors for codes from arbitrary algebraic curves. In their paper, 
Skorobogatov and ViaduJ also presented a modified algorithm, 
correcting more errors, but in general, not up to the designed 
minimum distance. Using profound results from algebraic 
geometry, Pellikaan [13] gave a decoding procedure which 
decodes up to [{d* - 1)/2J errors. However, his decoding 
procedure is very complex and is not completely effective. 
Recently, Justesen et at. [14] improved on their original 
decoding procedure in several ways and gave a new decoding 
procedure for codes from arbitrary regular plane curves, which 
can decode up to [(d' - g/2 - 1)/2J errot*. 

Manuscript received November 5. 1991; revised May 20. 1992. This work 
Mipponed in pari by the Office of Naval Research under Grant N000I4- 
»S I Hift?. This work was presented in pan at the 9th International Symposium 
tot Applied Algebra, Algebraic Algorithms, and Error-Correcting Codes, New 
Ottaito. LA, October 1991. 

Tte authors are with the Center for Advanced Computer Studies, University 
ot Southwestern Louisiana, Lafayette, LA 70504. 

IEEE Log Number 9203437. 



In ihis paper, we present a fairly simple decoding pro- 
cedure capable of decoding up to [(d* - 1}/2J errors. The 
improvement is obtained by using a form of majority scheme 
to find unknown syndromes in the well-known algorithm. The 
procedure can be implemented easily by hardware or software. 

The paper is organized as follows. In the next section, for 
easy reference, we include a fundamental iterative algorithm 
(FIA), which is very similar to the Gaussian elimination and 
can be used to easily derive the Berlekamp-Massey algorithm 
and the generalized Berlekamp-Massey algorithm [16]. Then 
we modify the FIA and give some related properties, which 
will be used in other sections. In Section III, a new decoding 
procedure for algebraic-geometric codes Ch(D, G) with G - 
mQ is presented. In order to easily understand this decoding 
procedure, one example is shown in Section IV. Finally, some 
conclusions are given in Section V. 

H. Fundamental Iterative Algorithm 
In this section, the fundamental iterative algorithm (FIA) 
[16] is modified. This modified algorithm is our main al- ' 
gorithm for decoding algebraic-geometric codes up to the 
designed minimum distance. To a certain extent it is similar 
to the Berlekamp-Massey algorithm, which is the main algo- 
rithm for decoding BCH codes up to the designed minimum 
distance. For easy reference, the FIA is described briefly in 
the following. This algorithm is for finding the smallest initial 
set of dependent columns in a matrix over any field F. That 
is, let 
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For each column j, let = ££o4*" IJ) **' 

where 4'"'"' ) = I. °e defined as the polynomial with the 
property that 
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for/iiSt-1, (2.2) 
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where o <ft) (:r) = 1 + a h ,ix + a h ,rx 2 + ••• + <ih„yx N and 
(/{*)]> expresses the coefficient of x j in f(x). Accordingly, 
the initial polynomial for column j is designated C ia '^(x) 
with C<°- l) (x) = 1 as the initial polynomial for the first 
column. Let 



(2.3) 

be the discrepancy at row i and column j. For some column 
j, if di.j a 0 for i = 1, 2, • • ■ , r- 1, then we have C^(x) = 
C<W)( X ) s ... a C( p - 1 -''>{.E)- Suppose d rJ # 0 and there 
is no «, where 1 < u < j and C M (x) = C (r - I,u >(i), 
then we define the final polynomial of column j, C u ^(x) = 
Ct r ~ l -fi{x), at the location of a rj - and place a cross "x" to 
indicate that the discrepancy is nonzero (we call it the final 
discrepancy of column j; the discrepancies at the locations of 
Oi,j for i < r are all zero), and then move on the next column. 
We then start examining the top element of the next column 
with CW<'-"(x) = C l '}{x) ~ C^-^ix). Otherwise, there 
is a C M {x) at a column u, where I < it < j such that 
&"Hx) = C^- l - n Hx). Then ^(i) can be obtained by 
the following equation 



C irJ) (x) = C^'^Hx) - 



i cM(x)xJ- v . 



The following can be easily seen. 

Lemma 1: There is an "x" in column j, if and only if it is 
linearly independent on its previous columns. 

Thus, beginning with the first column, we examine the 
elements in successive columns of matrix A, one by one, 
from the top row towards the bottom row. If the problem 
has a solution for A, C' (u ,, (a;) is a solution. From [16], it 
is easily seen that if </,,, = 0, then the subvector of the 
top i components of column j is a linear combination of 
the subvectors of the top * components of its previous j — I 
columns. Hereinafter, it is said that column j is a partial linear 
combination of its previous columns with the top i components. 
The linear combination coefficients arc the coefficients of 
CV<»{x). 

Thus, we have the following iterative algorithm, which is 
called the fundamental iterative algorithm (HA). Two storage 
tables, D and C, are set up for this algorithm. Table D stores 



the discrepancies rf r j, the final discrepancy of column j, while 
Table C stores the corresponding polynomials C (j) (x), the 
final polynomial of column j. 

Fundamental Iterative Algorithm: 

Step J) Empty Tables D and C, 1 -» j, 1 — r, 1 -» 

C(%). ; 
Step 2) Compute d rJ = [C^(x) ■ a^( x )}.. 
Step 3) IF d r j = 0, THEN 

a) IF r = M, then j - 1 -> /, C^(x) — C(x), 
stop; we have a solution. 

b) otherwise r + \—*r, and return to Step 2). 

Step 4) IF dr.j £ 0, THEN 

a) IF there exists a <f r ,u € D for some 1 < u < j, 
THEN 

dr.. 



C U) (x) - j± C^\x)x i - U ~» &i\x) 



(2.4) 



and return to Step 3). 
b) otherwise, d,.j is stored in D, C^(x) -» 
CU +i Hx) and C^(x) is stored in C, mark 
"x" at row r and column j. IF j < N, 
then j + I — j, 1 -+ r, and return to Step 
2). Otherwise stop, as this problem has no 
solution. 

Now let us consider another related problem. Let S' be such 
a matrix: (see matrix at bottom of page) where @ and # are 
ail unknown, and other elements are known. In each column 
j, if Si t j is @, then the values of s v j are known for u < t and 
s v j are all U for v > i; if there is no @ and Sij is the first 
#, then the elements above s. j are known and the elements 
below m,j are all #. Similarly, for each row i, if mj is an @, 
then the value of .•»,-,„ are all known for u < j, and $u, are all 
# for v > j; if there is no @ and » 4 j is the first #, then the left 
elements of ajj arc known and the right elements are all ft. 

Suppose column j is partial linear combination of its 
previous columns with the top i - 1 components and ay is an 
@, we want to know whether there is a unique value of @ at 
such that column j is a partial linear combination of its 
previous columns with the top i components. If there is, how is 
this unique value determined? The key step in decoding some 
linear error correcting codes can be reduced to this problem 
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mn j, while 

TU>(x), the 



T.) - C(X), 

o Step 2). 



(See next section): finding a partial linear combination with the 
top R components for a given integer R (if there is) and find all 
@, which are uniquely determined by the above requirement. 
In order to solve this problem, we modify the FLA as follows: 
2') The algorithm stops when <7 (R,j) (x) or CW{x) is 
obtained. 

2") Once C^ _1, ^(ar) is obtained, where s, 0 - is an @ or 
an #, then the algorithm defines the final polynomial 
of column j, C U) (x) = C7 (l " 1J >(i), and moves to the 
top element of the next column, namely, sij+i. 
Thus, after applying the modified FIA, the matrix D = 
(dij), called the discrepancy matrix, is obtained. The discrep- 
ancy matrix looks like the following: 
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From this matrix, it is easily to see the following. 

2a) If there is no "x" in column j, then column j is 
a partial linear combination of its previous columns 
with the top i - 1 components, where i is the smallest 
integer that an @ or # is at For the above matrix, 
i = 6, = 6; i - 4, j - 7; and so on. . 

2b) Suppose Sij is an @. This @ can be uniquely deter- 
mined such that column j is a partial linear combination 
of its previous columns with the top i components, if 
and only if there is no "x" at the locations of s,, u for 
1 < v < j. For example, @ at (4, 7) can be uniquely 



2c) If the value of @ at {i,j) is uniquely determined by 
(2b), then we have 



that is, 

Now we formally modify the FIA to solve this prob- 
lem. We need two more tables, E and F, which store these 
uniquely determined values of @ and their corresponding 
(7 (< - 3 ')(x), respectively. Let s< r )(x) = 1 + s r ,ix + 3 r ,2* 2 + • ■ • 

A Modified Fundamental Iterative Algorithm (MFfA): 



Step 1) Empty Tables D, C, E, and F, 

and 1 -» C7U>(z). 
Step 2) Compute d rJ = [C<'>(x)s< r >(x)] i . 
Step 3) IF d rJ = 0, THEN 



a) if r = R OR IF j = N and s r+liN is an @ or 
U, THEN STOP (according to (2')); 

b) if a r+ ],j is an @, r < R, and j < N, check if 
there is not s r+ i, u 6 D for 1 ^ u < j; when 
it is true, calculate -£fc»i c ft ) • Sr+ij-ftt 
this value and C^(x) are stored in E and 
F, respectively; then C U) (x) -» C^ vi >(x), 
jf + 1 -» j, 1 -» r, and return to Step 2) 
(according to (2") and 2a)-2c)); 

c) IF s r +ij- is an #, r < R and j < AT, then 
C<i)(a;) - C" +1 >(x), j + 1 - j, I - r, and 
return to Step 2) (according to 2'); 

d) otherwise, r + 1 -» r and return to Step 2). 

Step 4) IF d r , u ,<£ 0, THEN 

a) if there exists a dr iU e D for some 1< u < j, 
THEN 

<? W) (aO - ^C (u) (x)x'- U - C (i '(i) 

and return to Step 3); 

b) otherwise, d r j is stored in D, C^{x) is stored 
in C, mark "x" at the row r and the column 
j, C^(x) - C« +l >(x), J + 1 - j, 1 - r 
and return to Step 2). 

When the algorithm stops, a partial linear combination of 
column j with the top R components can be found or all 
values of @, which are uniquely determined by partial linear 
combinations, can be found. Obviously, the complexity of this 
• is 0(MN 2 ). 



III. A New Decoding Procedure 
In this section, we consider the general algebraic-geometric 
codes Crj(£), G) from X, where X is a nonsingular projective 

curve over F q , D is the divisor Pi H + P„, G = m • Q, 

and Q j4 P< for 1 < t < n. The linear code Cn(D,G) of 
length n over F, is the image of the linear map a* : Sl(G - D) 
F% defined by 

a'(u>) := (Resp 1 (w),Res ft (w),- -,Res P „(w)). 
From [15], the parameters of the code are 

d*=m-2ff+.2, k = m-g + l, 
t~[(d--l)/2\ = l(m + l)/2}- 9 . ■ 

We may assume that t > 0 or m > 2g. Recall that a number 
o, is nongap for Q if L(o,<?) 5^ Z-((oj - 1)Q). Then a function 
<f> 0i exists with pole order oi at Q. As is well known, the 
nongaps satisfy the following. 

Lemma 2: o 0 = 0, 

0<Oi <02 < ••• <o g -i < 2g, 
0{ = i + g, ' for i = g,g + !,•••, tn-g. 
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As a basis fi,h ; --- ,fm-g+i of the space L(G), we choose 
functions /,• = ^oj.u for t = 1,2, • • • ,rri + 1 — g. 

Let u, e, and c be a received word, an error vector, and a 
codeword, respectively. Then, we define the syndromes as 

Si := {«, = £ Uj&^i) - £ e <=< ** ) • < 3 - 1 ) 



where w is the number of errors. 
Now we define two-dimensional syndromes as 



It is known that if o,_i + Oj_i = o v ^ m, then 0 Oi _, • 
€ t(o^Q) C L(mQ). Thus, Sjj is a linear combination 
of .?o.5oit - , !S 0p and the coefficient of s 0p is not zero. 
Therefore, all Sij are known for Oj_i + o_,_i ^ m. 

On the other hand, Si.j may not be equal to s 0( _ 
because <j> Oi _ i ^ 0j _ l may not be equal to + 0j _, . However, 
there is a linear relation between them, that is, S,- j can be 
uniquely determined by s'k for 0 ^ k < + Oj_i. In 
this case, we say that s 0i ._ 1+0j _, and S,j are consistent and 
is a consistent term of s 0 j_ 1+0> _,. As a notation, we use 
5 o,_i+oj_, 10 express all consistent terms of s 0i _ l+OJ _ l . One 
example will be shown in Section IV. Hereinafter, "the value 
of ^Oi-i+oj-i" means "the value of +0j _,- and the values 
of its possible consistent terms"; "the number of a 0i ,, +0j _, 
in matrix S*" means "the number of it and its possible 
consistent terms in S*"; "to substitute @ for a Oi _ l+0j _ l " 
means "to substitute @ for it and its possible consistent terms"; 
"the number of s 0i-l + 0 j.i with the same value" means "the 
number of it and its possible consistent terms, whose values 
satisfy some liner relations"; and so on. 

We construct a matrix, S =(Si./)( m _ 9+ i )x(m _ 9+1) . from 
Lemma 2: 



5,,,-,., S m - 9 . 

.Sm-9+l.l Sm-j+ 



Sl.3 
St.* 
I 

5m-,,: 



••• Si.m-s $l t m-,+l 
■•• Si.m-a Sj.m-j+l 
I • I 

— S m -,. m - 9 S m - s . m . s+l 

•••Sm-f+l.m-3 Sm-»+l.m-j+l 



1 

1 0„,(P ta ) 

I I 

1 <*»,(PtJ 



I • 



<><.„,_„_ .(P*,) <j„ m _,(P tl ) 

Oo M _,_,(Pt s ) <v„_,(Pt a ) 

I I 

<*»„_,-. <P*J *.„_,( Pi.). 



Suppose the value of in S are all known, we have the 
following theorem. 

Theorem 1: If column j of S is a partial linear combination 
of its previous columns with the top [(m + 1)/2J components, 
and if the coefficients are a,- for 1 < t < j — 1, then column 
j is linearly dependent on its previous columns and the error 
rational points are the roots of 



(3.2) 



/> -£a,-./ i = 0. 



Proof: See [12, Theorem 1J. 



(3.3) 



Unfortunately, the values of s m +i , • • • , s m+g are unknown, 
that is, the values of Sij for a-i + oj-x = m + 
l,m + 2, • • • ,m + g ate unknown, (3.3) may not be found 
from the matrix S. Thus, the key problem in decoding 
algebraic-geometric codes is finding of the real values of s, 
for t = rn + 1. m + 2, ■ • • , ro + q and (3.3). This problem 
can be solved by the MFIA developed in Section II. In order 
to find the values of s,- for j = m + 1, m 4- 2, ■ • • , m + g, 
our decoding procedure is to find them iteratively, i.e., first 
find s m+ i, then s, n+2 . and so on, till s m+g is found. In the 
following, we describe how to find s m+w if we know s, for 
i = O.Oi,-- -,rrv -,771 + u; - 1, where 1 $ w ^ g. 

Suppose s m+ \..--- ,s m +w-\ are found and s m + w is 
still unknown, where 1 < w < g. Now we want to find 
3 m +w Let 5* be the rewritten S, where s p is substituted for 
Si,j{oi-i + Oj_i = p) for p < m + w - 1, @ is substituted 
for S u ,t,(Oii-i + o«-i = m + w), i.e., for 3 m+tu , and # is 
substituted for $k,h(ok-i + o/,_j > m + w), i.e., for s, with 
q > m + w. Obviously, S' satisfies the requirements of S' 
in Section II. We will use the MFIA to find the value of @, 
i.e., the value of s m + w . 

Before describing how to find the value of $ m+w , we first 
calculate the number of @, i.e., the number of a m+w in 5*. 
The next lemma is useful. 

Lemma 3: Let A w be the number of nongaps in [1, w - 1) 
for 1 < to < g, then 



ii.<l(«-l)/2J. 



(3.4) 



From (3.1) and (3.2), Sij = s 0i _, and S.-.i = «„,._,, and S can 
be decomposed into a product of three matrices as follows: 



Proof: Obviously, it is true for w = 1. For 2 < w ^ g, 
e consider the two cases 

Case 1) If w is a gap, then we consider the two subcases. 

a) w is odd. If s € [l,tu - 1], then w - s € 
[1. w — 1 j and s t^w — s, s and «; — s are not 
both nongaps. Thus, (3.4) is true. 

b) w is even, w/2 € [l,w — 1] can not be a 
nongap. If s € [1, w - Ij and a ^ w/2 then 
tu — 3 6 [1, w — 1], a 5^ w — a, s and u> — a 
are not both nongaps. Thus, A w ^(w- 2)/2 
and (3.4) is true. 

Case 2) If w is a nongap, then from Lemma 2 and 2 < w ^ 
g, we can assume that w, w+ 1, ■•■,w + p- 1 are 
nongaps and to+p is a gap, where 1 < p < g. From 
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j the previous proof, A m+P ^ [(w + p - 1)/2J. 

On the other hand, A wrP = A,„ + p. we have 
I A m + p [(w + p - 1)/2J and (3.4) is true. □ 

J Thus, for the number of @, we have the following theorem. 
) Theorem 2: The number of @ in S* is at least m - 2g, thai 
j is, at least d* - 2. 

Proof; For each l<j<m-g+l, from (3.1) there is 
an @ in column j, if and only if 

.m + w-Oj-i £ {o 0 ,Oi,---,o m _ a }. (3.5) 
Equivalently, there is no @ in column j, if and only if 

m + w-Oj-i g {o 0 ,Oi,-- ,o m _ 9 }. (3.6) 

Let us calculate the number of j, which satisfy (3.6). We 
consider the following two cases 

a) If m + w - Oj_i > m, i.e., Oj_i < w, then (3.6) is 
satisfied. From Lemma 3 the number of such j is l+A w 
0 = 1 satisfies this condition). 

b) If m + w - Oj_i ^ m, i.e., Oj_j ^ w, then the num- 
ber of to + w - Oj-\ being gaps is B w , which is the 
number of the gaps in [w, 2g - 1). 

. From this discussion, the number of columns, in which there 
are no @, is 1 + A w + B w . Since the number of nongaps in 
[Lw-l] isA w , the number of gaps in [l,u/-l] is u>— 1 — v4„,. 
From Lemma 2, we have 

w-l-A w + B w = g. (3.7) 
Thus, from Lemma 3, we have 

1 + A v + B w = 1 + A w + g - w + 1 + A v < 1 + g, 



s at least (m - g + 1) - (g + I) = 



namely, the number of <S 
m - 2g. U 
Now we are going to explain how to use the MFIA to 
find the value of s m+ „,. Some useful facts are introduced as 
follows. 

Lemma 4: If (A) there is a unique value of this @ such that 
column j of S* is a partial linear combination of its previous 
columns with the top i components and (B) column j of S is 
linearly dependent on its previous columns, then the unique 
value must be equal to s m+V! . 

Conversely, from Lemma 1 and Lemma 4, we have the 
following. 

Lemma 5: If (A) is true and (C) the unique value is not 
equal to s m+ul , then column j of S is linearly independent 
from the previous columns, and when the FIA is applied to 
S, there is an "x" at 

Fdr convenience, if (A) is true, then the unique value is 
called as a candidate value for s m+w , wr simply, candidate. 
Thus, in order to determine s m+w , our first objective is to find 
all candidates. In Section II, we have developrd the MFIA to 
find all candidates. If there is an "x" at we say this "x" 
affects @ in row i and @ in column j. When the MFIA is 
applied to S*, we can see the following properties. 

a) If there is an "x" in column j, then column j of S is 
linearly independent on its previous columns. 



b) If an @ is at (i, j), then it can be uniquely determined by 
a partial linear combination of its previous columns with 
the top i components, if and only if no "x" affects it. 

c) An "x" at (1, 1) does not affect any @. An "x" in the 
first row and not in the first column affects at most one 
@. An "x" in the first column and not in the first row 
affects at most one @. An "x" neither in the first row 
nor in the first column affects at most two @. 

Now we have the following theorem. 

Theorem 3: Applying the MFIA to matrix S*, there is at 
least one @, which can be uniquely determined by a partial 
linear combination, namely, there is at least one candidate 
value for s m+w calculated by Step 3b). 

Proof: Let t = [(a" - 1)/2J. Since v < t, there are at 
most v linearly independent columns in X as well as in S. 
Suppose there are \i 14 x" in the discrepancy matrix D after 
applying the MFIA to S\ From Property a), /j ^ v sS t. We 
consider the following three cases. 

1) If an "x" is at (1, 1), then it does not affect any @, and 
the other }i - 1 "x" affect at most 2(p - 1) @ from 
Property c). Thus, at most 2/x - 2 @ are affected. 

2) If an "x" in the first column is not in the first row, then 
it affects at most one @ from Property c). But, there 
must be another "x" in the first row; this "x" affects at 
most one @, too (if the received vector has v errors for 
1 $ v ^ t, then the known syndromes are not all zero 
and there must be one "x" in the first row and in the 
first column, respectively). Thus, the other fj. - 2 "x" 
can affect at most 2(/i - 2). Hence, at most total 2p. - 2 
@ are affected. 

3) If p. = 1, from the discussion in 2), this "x" must be 
at (1, 1). Thus, there is no @ affected from Property c), 
that is, 2/i - 2 (in this case, 2/< - 2 = 0) @ are affected. 

From Theorem 2, the number of @ is at least d r -2, at least 
d" - 2 - (2/i - 2) ^ 1 @ is not affected and it is a candidate 
from Property b). □ 

Generally speaking, after applying the MFIA to S* there 
are often more than one candidate and they may be not all 
correct. Fortunately, s m + w can be easily determined from all 
candidates by the following theorem. 

Theorem 4: After applying the MFIA to 5*, the number of 
correct candidates is greater than the number of the incorrect 
candidates. 

Proof: Suppose that there are n "x" in the discrepancy 
matrix D after applying the MFIA to S* . From the proof of 
Theorem 3, we can obtain at least d* - 2 - (2/i - 2) candidates. 
If a candidate is equal to s m+m , it is called a correct candidate, 
otherwise it is called an incorrect candidate. Now we are going 
to calculate the number of incorrect candidates. 

From Property b), the columns of S, which correspond to 
these columns of S' containing "x", are linearly independent 
on their previous columns. The number of these columns is 
/t. Since the total number of linearly independent columns of 
S is v or less (because the rank of X is v or less), and from 
Lemma 5, the number of the incorrect candidates is at most 
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v ~ ix. Therefore, at least d* - 2 - 2(/i - I) - (v - ft) = be a basis of L(G). Then, 
d' - v — /t candidates are correct, and the number of the 
correct candidates is greater than the number of the incorrect 
candidates (d* - v - ft > v - n). The proof is completed. □ 



Once s m+U! is found, S* is modified: substituting the correct 
value of s m+w for @ and substituting @ for s m+w +i. Since 
the number of .? m+u ,. + i, that is, of new @, is d' ~ 2 or more, 
if (3.3) is not found, then by the same reason and in the same 
way, » m+a . +1 can be found. 

However, it is easily seen that to apply the MFIA to the 
new S* is equivalent to modifying the discrepancy matrix 
D as follows: substitute the value of a m + w for all @ being 
not candidates and eliminate the discrepancy at this place if 
possible, that is, there is no "x" on its column and there is one 
"x" on its row; substitute 0 for @ being a correct candidate 
and "x" for @ being a incorrect candidate (by Lemma 4 and 
Lemma 5); and place @ at (ij), at which s m+w+ i was. The 
complexity of modifying discrepancy matrix D is 0((m — 
g + l) 2 ). Thus, in our decoding procedure the construction of 
5* has complexity 0((m -g+l) 2 n). The complexity of the 
algorithm for determining value of s m+ i is 0((m-g+l) 3 ). To 
find « m+ 2, a m+ 3, • • • , s m+g at most g - 1 modifications of D 
are required (for some error patterns, (3.3) can be early found), 
which has complexity 0((g - L)(m - g + l) 2 ). Therefore, to 
find all Si [for i = m+ 1, ■ • • , m +g, and (3.3), the complexity 
is 0((m -g + l) 2 n). 

Our decoding procedure can be outlined as follows. 

1) Calculate the syndromes Si from the received vector for 
i = oq,oi,--- ,o m - a+ i. 

2) Determine s,- for » = m + 1, • • • , ro + g and (3.3). 

3) Find the roots of (3.3) among the rational points. 

4) Solve a system of linear equations and obtain the error 
locations and error magnitudes at the same time. 

Obviously, this decoding procedure is a generalization of 
Peterson's decoding procedure for the BCH codes. It is easily 
seen that the complexity of this decoding procedure is 0(n 3 ) 
where m ~ n. 



IV. An Example 
In this section through an example, we show this simple 
decoding procedure. Let us consider a plane curve defined 
by the equation f(X,Y,Z) - X 6 + Y 4 Z + YZ 4 =0 
over GF(2 4 ). This curve is called the Hermitean curve. With 
Q = (0 : 1 : 0), x = X/Z, has pole order 4 and y = Y/Z 
has pole order 5. From [15] we know that Hermitean curve 
has the genus g = 6 and 65 rational points in GF(2 4 ). Define 
the code with length 64, i.e., D is all finite rational points, 
and G = 23Q. The linear code C n (D,G) of length 64 over 
GF(2 4 ) is the image of the linear map a' : Sl{G - D) -» 
GF{2 4 ) M defined by 

a*(w) := (Res Pl (u/),R<s^(u;),-.-,Resp 6 >)). 

From Riemann-Roch theorem, the designed minimum dis- 
tance d* = deg(G) - 2g + 2 = 13. Any six or fewer errors 
can be corrected. Since deg(G)-#+l = 18, let /i,/a,-*-,/i8 
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h{Px) fi(Pi) 

| | ... | I (4 .1) 



is a parity check matrix for Cn(D,G). 
A basis of L{G) is 

/i = l;(0) / 2 = x;(4) / 3 = y;(5) 
/ 4 =* 2 ;(8) fs = xy;(9) / 6 = y 2 ;(io) 
/ 7 = r»;(12) / 8 =s 2 i/;(13) / 9 =xj/ 2 ;(14) 
/io = y 3 i (15) fn - i 4 ; (16) /i2 = x 3 y; (17) 
/i3=*V;(18) / 14 = x2/ 3 ;(19) / 15 =y 4 ;(20) 
/ 16 =x 4 j/;(21) /i7=iV;(22) /is = xV;(23) 

(4.2) 

where the number in parenthesis indicates the pole order of 
the corresponding function with Q. 

For convenience, let 4>a = /i. 04 = h, 05 = h, 4>& — Si, 
<h = h, ho = /e. and fa = /<_ 5 for i = 12, 13, • • ■ , 23. 
Thus, 



ordQ(&) = 



= 0,4,5,8,9, 10,12, 13,- -,23. 



For each i, there may be many 0 such that ordQ(0<) = ~i. 
But it is sure that these 0 must be linear combinations of 
#o>^3, &»"■.&. For example, ordcj(i 5 ) = -20. Since 
X s + Y 4 Z + YZ 4 = 0, we have 

a; 5 = y 4 + y, 
x 6 = xy 4 + xy, 

and so on. 

Let 020 = y* and 0 2o = x 5 , 0 2 4 = xy 4 , 0 24 = x 6 , and so 
on. Then, from the previous equations, we have 

020 = 020 + 05. 



024 = 024 + 09- 



(4.3) 



Let w = (tti,ti2,-"i"64) be a received word, e = 
(ei,e2,--,e64) be an error vector, and c= (ci,C2,--,C64) 
be a codeword. Thus, we have u = e + c Then, we define 



:= («,a*0i) 

=f;«,-0 i (p j ), 



for 1 = 0,4,5,8,9, 10, 12, 13, • 



(4.4) 



64 

s'i := («,a*0j) = E ^W)." for f = 20 
= Ee^i). (4-5) 

64 64 

ai = E c i^.( p i) a nd *< = E e J0<( p y) f( > r «' > 24 > 
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(4.1) 



(5) 

;(io) 
2 ;(H) 
3 y;(ir) 
*; (20) 
V;(23) 

(4.2) 



where 0,- and $ have order -i with Q. Certainly, s< and $' f 
for i ^ 24 are unknown. From (4.3)-(4.5), we know that 



Suppose there are v errors e*„ at Pk^ = (At,. : : Zjt„ ) for 
H = l,2,---,v, where v ^ 6. 
Consider a received vector 



From (4.2) and (4.4), we have 



u = (ft 12 . a\ a 7 : a 8 , ft 9 , « 9 , 0, • • • , 0, 0, 0, 0, • • • , 0, 
0,0,0.0,0), 

that is, mi = ft 12 . «2 = a 4 , «3 = a 7 , ttj = a 8 , us = a 9 , 
w B = « 9 , and other components are 0. Let Pi = (1 : 1 : a), 
P 2 = (1 : 1 : a 2 ), P 3 = (1 = 1 : « 4 ). P< = (1 : 1 : « 8 ), 
P 5 = (0 : 0 : 1), Pe = (1 : « : 1). 



s,8 = « 8 sio = a 2 

S 2 2 = ft 5 3 2 3 = ft 1 ' 



.916 = 1 r 
.120 = ft 1 



Si3 = a 11 

S17 = ft 6 



In order to find 324 (and s' 24 ), we construct S* as shown in 
the first matrix at bottom of page where @ expresses s 2 4 or .s 24 
and # expresses other unknown syndromes. In this example, 
S* is shown in the second matrix at bottom of page. 

Applying the MFIA to 5*. we have the following discrep- 
ancy matrix, where the candidate in column j is denoted by @j 
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and "x" in column j is denoted by Xj, and their values arc in 
the first matrix at bottom of page where xi = a, x 2 = 1, and 
x 3 a a 5 ; «U = «- = ®u = ®s = «'«o = «*a = ^10 = 
a 8 . From Theorem 4, .-^ = 8S3 = ® 6 = #9 = ««io = ft* 
are correct and «^ = = <$n = rt H are not correct. $04 
should be « 3 . 

Now we modify the discrepancy matrix as follows. Substi- 
tute n H (or ft 3 ) for all @ being not candidates and eliminate the 
discrepancies at these places if it is possible (the complexity is 
at most 0(t?), then substitute @ for .125 (or .i 25 ). Then, place 
0 at (5,10), (6,9), (9,6), and (10,5), respectively (the candidate 
values for k-za at these locations were correct) and place x^ at 
(11,4), x 7 at (7,7), and Xu at (4, 11), (the candidate values 
for s 24 arc not correct). The new matrix is equivalent to that the 
MFI A has been applied to new 5*. where the values of .s 2 j and 



S24 are known and the values of 035 and s 25 are unknown. So 
to find *25 and d' 25 , that is, calculate the values of candidates of 
all .s'25 and eliminates of discrepancies at some places, where 
»2A were. After this discrepancy matrix is modified, we have 
the second matrix at bottom of page where Xi = et ( x 2 = 1, 

X 3 = ft 5 , X A = X 7 = X], = ft 12 ; ,?2h as ©e = ®io = 1. 

By the same way we can find that * 2 e = a, .S27 = a 4 , 
82a = 1, and .s 2 9 = « 8 . Thus, we can find two solutions of 
(3.3). They are linear dependent relations of column 5 with • 
the first 4 columns and column 6 with the first 4 columns, 
respectively. We have 



' xy + X? + y + X as 0. 
y 7 + x 2 + + <t*y = 0. 
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known. So 
ndidates of 
ces, where 
1, we have 
», x 2 = I, 
u = 1. 
* w = u 4 , 
alutions of 
mn 5 with 
1 columns, 



The two equations have six common roots, the five points on 
the line y — x: 



^ = {1 
P S = (0 



:l:o»), 
:l:a 8 ), 



and the point of intersection of the two lines x + 1 = 0 and 

y + x + « 4 : 

F 6 = (1 : a : 1). 

in order to find the real error locations and the error 
magnitudes, we should solve reduced parity check equations 
in the six unknowns e«, Ci, e 3 , e 4 , e 5 , and ee. Solving these 
linear equations, we have 



Thus, there are six errors. They are at P\, Pi, Pz, Pi, P&, 
and P$ and the error magnitudes are r* 12 , a 4 , a 7 , a 8 , a 9 , and 
n; 9 , respectively. 

. V. Conclusion 
In this paper, we have derived a very simple decoding 
procedure for decoding algebraic-geometric codes Cq(D,G) 
with G = m • Q up to [(d* - 1)/2J errors. The computa- 
tion complexity of this decoding procedure is 0(n 3 ). This 
decoding procedure is a generalization of Peterson's decoding 
procedure. It should be noted that the decoding procedure 
is applicable to decoding some cyclic codes up to the van 
Lint-WHson's bound and decoding some cyclic codes beyond 
actual minimum distance. 
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